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Abstract. We present an updated analysis of the constrained minimal supersymmetric 
standard model with /i > supplemented by an 'asymptotic' Yukawa coupling quasi-unification 
condition, which allows an acceptable b-quark mass. Imposing constraints from the cold dark 
matter abundance in the universe, B physics, the muon anomalous magnetic moment, and 
the mass mj, of the lightest neutral CP-even Higgs boson, we find that the lightest neutralino 
cannot act as a cold dark matter candidate. This is mainly because the upper bound on the 
lightest neutralino relic abundance from cold dark matter considerations, despite the fact that 
this abundance is drastically reduced by neutralino-stau coannihilations, is incompatible with 
the recent data on the branching ratio of B s — > fj, + . Allowing for a different particle, such as 
the axino or the gravitino, to be the lightest supersymmetric particle and, thus, constitute the 
cold dark matter in the universe, we find that the predicted m^'s in our model favor the range 
(119 - 126) GeV. 



1. Introduction 

The well-known constrained minimal supersymmetric standard model (CMSSM) [H El El H]; 
which is a highly predictive version of the minimal supersymmetric standard model (MSSM) 
based on universal boundary conditions for the soft supersymmetry (SUSY) breaking parameters, 
can be further restricted by being embedded in a SUSY grand unified theory (GUT) with a gauge 
group containing SU{A) C and SU(2)r. This can lead [5] to 'asymptotic' Yukawa unification 
(YU) [6j, i.e. the exact unification of the third generation Yukawa coupling constants ht, h^, 
and h T of the top quark, the bottom quark, and the tau lepton, respectively, at the SUSY GUT 
scale Mqut- The simplest GUT gauge group which contains both SU(A) C and SU(2)r is the 
Pati-Salam (PS) group G PS = SU(4) C x SU(2) L x SU(2) R E] - for YU within 50(10), see 
Refs. [91 [TO]. 

However, given the experimental values of the top-quark and tau-lepton masses (which, 
combined with YU, naturally restrict tan j3 ~ 50), the CMSSM supplemented by the assumption 
of YU yields unacceptable values of the 6-quark mass mj, for both signs of the parameter [i. This 
is due to the presence of sizable SUSY corrections [11] to mj, (about 20%), which arise [TT| [12] 
from sbottom-gluino (mainly) and stop-chargino loops and have the sign of /i - with the standard 
sign convention of Ref. [13]. The predicted tree-level rrib(Mz), which turns out to be close to 
the upper edge of its 95% confidence level (c.l.) experimental range, receives, for fi > [fi < 0], 
large positive [negative] corrections which drive it well above [a little below] the allowed range. 



Consequently, for both signs of fi, YU leads to an unacceptable mb(Mz) with the \i < case 
being much less disfavored. 

The usual strategy to resolve this discrepancy is the introduction of several kinds of 
nonuniversalities in the scalar [91 [TO] and/or gaugino |14[ [T5] sector of MSSM with an 
approximate preservation of YU. On the contrary, in Ref. |16j . concrete SUSY GUT models 
based on the PS gauge group are constructed which naturally yield a moderate deviation from 
exact YU and, thus, can allow acceptable values of the 6-quark mass for both signs of [i within the 
CMSSM. In particular, the Higgs sector of the simplest PS model [5] is extended so that the 
electroweak Higgs fields are not exclusively contained in a SU(2)l X SU(2)r bidoublet superfield, 
but receive subdominant contributions from other representations too. As a consequence, a 
moderate violation of YU is naturally obtained, which can allow an acceptable fo-quark mass 
even with universal boundary conditions. It is also remarkable that the resulting extended SUSY 
PS models support new successful versions [T7] of hybrid inflation based solely on renormalizable 
superpotential terms. 

These models provide us with a set of 'asymptotic' Yukawa quasi-unification conditions, which 
replace exact YU and can be applicable [HI [181 02] for both signs of the MSSM parameter fi. 
We focus here on the (i > case since \i < is strongly disfavored by the constraint arising 
from the deviation 5a^ of the measured value of the muon anomalous magnetic moment 
from its predicted value a^ M in the standard model (SM). Indeed, \i < is defended only at 3a 
by the calculation of a^ M based on the r-decay data which is presented in Ref. |20| . whereas 
there is a stronger and stronger tendency at present to prefer the e + e~-annihilation data for the 
calculation of a^ M , which favor the \i > regime. Moreover, in Ref. [2T], it was claimed that, 
after some improvements, the r-based result shifts considerably towards the e + e~-based one. 

Let us recall that, in this case, the suitable 'asymptotic' Yukawa quasi-unification condition 
applied [MUHICE] is 

h t :h b :h T = \l + c\:\l-c\:\l + 3c|. (1) 

This relation depends on a single parameter c, which is taken, for simplicity, to be real and 
lying in the range < c < 1. With fixed masses for the fermions of the third generation, we 
can determine the parameters c and tan (3 so that Eq. (fTJ) is satisfied. In contrast to the original 
version of the CMSSM [2j [31 0], therefore, tan (3 is not a free parameter, but can be restricted, 
within our set-up, via Eq. ([1]) to relatively large values. The remaining free parameters of our 
model are the universal soft SUSY breaking parameters defined at Mgut> i>e. 

Mi/2, rn , and A , (2) 

where the symbols above denote the common gaugino mass, scalar mass, and trilinear scalar 
coupling constant, respectively. These parameters can be restricted by employing a number 
of experimental and cosmological requirements as in Refs. [16} \TE[ [TO] and most recently in 
Refs. [22^ 123] . In this talk, we review the results of Ref. [22J implementing the following 
improvements: 

• We do not take into account the upper bound on rriy implied by the lower bound on Sa^ 
from the r-based calculation of Ref. [20] raising, thereby, the upper bound on my from the 
muon anomalous magnetic moment - see Sec. [2j 

• We employ the recently released data on the branching ratio of B s — > fi + fj>~ [24J and the 
mass rrih of the lightest CP-even Higgs boson [25]. As a consequence, our predictions in 
Ref. [22j for rrih and the role of the lightest neutralino as cold dark matter (CDM) particle 
have been significantly altered. 

All the cosmological and phenomenological requirements which we considered in our 
investigation are exhibited in detail in Sec. [2j Restrictions on the parameter space of our model 
are derived in Sec. [3] and our conclusions are summarized in Sec. [H 



2. Cosmological and Phenomenological Constraints 

In our investigation, we integrate the two-loop renormalization group equations for the gauge 
and Yukawa coupling constants and the one-loop ones for the soft SUSY breaking parameters 
between Mqtjt and a common SUSY threshold Msusy — ( rn t 1 m t 2 ) 1 ^ 2 (*i f 2 are * ne S ^°P mass 
eigenstates) determined in consistency with the SUSY spectrum. At Msusy , we impose radiative 
electroweak symmetry breaking, evaluate the SUSY spectrum employing the publicly available 
calculator SOFTSUSY [26j, and incorporate the SUSY corrections to the b and r mass [12]. The 
corrections to the r-lepton mass m T (almost 4%) lead [161118] to a small decrease of tan /3. From 
-Msusy to Mz , the running of gauge and Yukawa coupling constants is continued using the SM 
renormalization group equations. 

The parameter space of our model can be restricted by using a number of phenomenological 
and cosmological constraints. We calculate them using the latest version of the publicly available 
code micrOMEGAs [27J. We now briefly discuss these requirements - for similar recent analyses, 
see Ref. [28] for the CMSSM or Refs. QH [29] for the MSSM with YU. 

SM Fermion Masses. The masses of the fermions of the third generation play a crucial role in 
the determination of the evolution of the Yukawa coupling constants. For the 6-quark mass, we 
adopt as an input parameter in our analysis the MS 6-quark mass, which at la is |30] 

m h (m b ) m = 4.19±8;S| GeV. (3) 

This range is evolved up to Mz using the central value a s (Mz) = 0.1184 [30J of the strong fine 
structure constant at Mz and then converted to the DR scheme in accordance with the analysis 
of Ref. [31] • We obtain, at 95% c.L, 

2.745 < m b (M z )/GeV < 3.13 (4) 

with the central value being m&(Mz) = 2.84 GeV. For the top-quark mass, we use the central 
pole mass (M t ) as an input parameter [32] : 

M t = 173 GeV =>■ m t (m t ) = 164.6 GeV (5) 

with nit(mt) being the running mass of the t quark. We also take the central value m T {Mz) = 
1.748 GeV [31] of the DR tau-lepton mass at M z . 

Cold Dark Matter Considerations. According to the WMAP results [33], the 95% c.l. range 
for the CDM abundance is 

Q CDM h 2 = 0.1126 ±0.0072. (6) 

In the context of the CMSSM, the lightest neutralino x can be the lightest supersymmetric 
particle (LSP) and, thus, naturally arises as a CDM candidate. In this case, the requirement 
that the x relic abundance Q^h 2 does not exceed the 95% c.l. upper bound from Eq. ([6]), i.e. 

n^h 2 < 0.12, (7) 

strongly restricts the parameter space of the CMSSM. This is because Q^h 2 increases, in general, 
with and so an upper bound on can be derived from Eq. ([7]). The calculation of Q,^h 2 
in micrOMEGAs includes accurately thermally averaged exact tree- level cross sections of all the 
possible (co) annihilation processes [31 [M], treats poles [H US [35] properly, and uses one- loop 
QCD and SUSY QCD corrections [HI [161 ES] to the Higgs decay widths and couplings to 
fermions. It should, though, be noted that the restrictions induced by Eq. ([7]) can be evaded if 
we adopt one (or a combination) of the following scenarios: 



• The cosmological evolution before Big Bang nucleosynthesis (BBN) deviates from the 
standard one \37\ 138] , Since x within the CMSSM is essentially a pure bino, the scenario 
which fits better this case is the low reheat temperature scenario with the decoupling of 
X occurring before reheating. This scenario, however, is disfavored since it requires a very 
low reheat temperature ~ (1 — 5) GeV. We will, thus, assume that the decoupling of the 
neutralino from the cosmic fluid occurs during the conventional radiation dominated era. 

• The lightest neutralino is not the LSP and, thus, the relic density of another SUSY particle 
[39] . which is the LSP, is to account for f^cDM^ 2 - This particle could be the gravitino [30] 
or the axino [41[ 142] . The case of gravitino is tightly restricted in the CMSSM due to the 
BBN constraints imposed during the decay of the lightest ordinary super 'symmetric particle 
(LOSP) to it. On the other hand, axino CDM [42] is, in general, possible once its mass 
and the reheat temperature are chosen appropriately. In such a case, x may play the role 
of the LOSP and can contribute to the non-thermal production of the LSP. In particular, 
its contribution to the relic density of the LSP is equal to Q^h 2 times the ratio of the LSP 
mass to my. and, thus, f^/i 2 's exceeding the bound in Eq. © can be perfectly acceptable. 

The Branching Ratio BR (6 —¥ sj) of b — >■ sj. The most recent experimental world average 
for BR(6 -)• sj) is known @3] to be (3.52 ± 0.23 ± 0.09) x 10~ 4 and its updated SM prediction 
is (3.15 ±0.23) x 10 -4 |44j . Combining in quadrature the experimental and theoretical errors 
involved, we obtain the following constraints on this branching ratio at 95% c.l.: 

2.84 x 10" 4 < BR(6^s 7 ) < 4.2 x 10" 4 . (8) 

The computation of BR (6 — > S7) in the micrOMEGAs package presented in Ref. [36] includes |45j 
next-to-leading order (NLO) QCD corrections to the charged Higgs boson (H^) contribution, 
the tan/3 enhanced contributions, as well as resummed NLO SUSY QCD corrections. The 
contribution interferes constructively with the SM contribution, whereas the SUSY contribution 
interferes destructively with the other two contributions for [i > 0. The SM plus the and 
SUSY contributions initially increases with and yields a lower bound on from the lower 
bound in Eq. ([5]). (For higher values of m^, it starts mildly decreasing.) 

The Branching Ratio BK(B S — > fj> + fi~) of B s — > ^ . The rare decay B s — > /x + occurs via 
Z penguin and box diagrams in the SM and, thus, its branching ratio is highly suppressed. The 
SUSY contribution, though, originating [46 1 137] from neutral Higgs bosons in chargino-, 
and W^-mediated penguins behaves as tan 6 /3/m\ (itia is the mass of the CP-odd Higgs boson 
^4) and hence is particularly important for large tan /3's, especially after the new more stringent 
95% c.l. upper bound 

BR (B s < 1-08 x 10~ 8 (9) 

recently reported by CMS and LHCb [23]. This new bound significantly reduces the previous 
bound [35], which we had adopted in Ref. [22]. The bound in Eq. ([9]) implies a lower bound on 
niy since BR(i? s — > /i + /x~) decreases as ttt-lsp increases. 

The Branching Ratio BR {B u — > tu) of B u — >■ tu. The purely leptonic decay B u — > tu proceeds 
via W^ 1 - and ff^-mediated annihilation processes. The SUSY contribution, contrary to the SM 
one, is not helicity suppressed and depends on the mass m H ± of the charged Higgs boson since 
it behaves (32 [39] as tan 4 /?/™ 4 ^. The ratio R(B U -> tu) of the CMSSM to the SM branching 
ratio of B u — > tu increases with and approaches unity. It is to be consistent with the 
following 95% c.l. range [33] : 

0.52 < R (B u -> tu) < 2.04 . (10) 
A lower bound on can be derived from the lower bound in this inequality. 
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Figure 1. Summary of the conventions 
adopted in Figs. [2] and [3] for the various 
restrictions on the model parameters. 



Muon Anomalous Magnetic Moment. The quantity 5a^, which is defined in Sec. [H can be 
attributed to SUSY contributions arising from chargino-sneutrino and neutralino-smuon loops. 
The relevant calculation is based on the formulas of Ref. [50J. The absolute value of the result 
decreases as increases and its sign is positive for /j, > 0. On the other hand, the calculation of 
not yet completely stabilized mainly because of the ambiguities in the calculation of the 
hadronic vacuum-polarization contribution. According to the evaluation of this contribution 
in Ref. [2D], there is still a discrepancy between the findings based on the e + e~-annihilation 
data and the ones based on the r-decay data - however, in Ref. [2T], it was claimed that this 
discrepancy can be considerably ameliorated. Taking into account the more reliable calculation 
based on the e + e~ data and the experimental measurements [JT] of a«, we obtain the following 
95% c.l. range - cf. Ref. [52] : 

12.7 x 1(T 10 < <5a M < 44.7 x 1(T 10 . (11) 

A lower [upper] bound on can be derived from the upper [lower] bound in Eq. (|11|) . As 
it turns out, only the upper bound on is relevant in our case. Taking into account the 
aforementioned computational instabilities and the common practice [2H] , we consider this bound 
only as an optional constraint. 

Collider Bounds. For our analysis, the only relevant collider bound is the 95% c.l. LEP bound 
|53] on the lightest CP-even neutral Higgs boson mass 

m h > 114.4 GeV, (12) 

which gives a lower bound on m^. However, we should keep in mind that recent data from 
ATLAS and CMS [25] provide a 99% c.l. upper bound < 128 GeV and a hint in favor of the 
range (125 ± 1) GeV. Allowing for a theoretical error of ±1.5 GeV and adding in quadrature 
the experimental and theoretical uncertainties, we construct the la range of interest [54]: 

123.2 < m h /GeY < 126.8. (13) 

The calculation of m;, in the package S0FTSUSY [26] includes the full one-loop SUSY corrections 
and some zero- momentum two- loop corrections [55] • The results are well tested [56] against 
other spectrum calculators. 



3. Restrictions on the SUSY Parameters 

Imposing the requirements described above, we can delineate the allowed parameter space of 
our model. Throughout our investigation, we consider the central values for the SM parameters 



M t , mb(Mz), m T (Mz), and a s (Mz)- We adopt the following conventions for the various lines 
and regions in the relevant figures (Figs. [2] and [3]) - see Fig. [TJ 

• On the solid black line, Eq. ([7]) is saturated. 

• The horizontally hatched region is allowed by Eq. ([7]). 

• In the light gray region, the lightest stau r 2 is lighter than x- 

• The dark gray region is excluded by the lower bound in Eq. ([8]). 

• The gray region is excluded by Eq. ([9]). 

• The yellow region is excluded by the lower bound in Eq. ([TO]) . 

• The vertically hatched region is favored by the lower bound in Eq. 

• The red region is excluded by Eq. (fT2j) , 

Note that the upper bounds in Eqs. ([8]), (fTUj) . and (fTTj) do not restrict the parameters of our 
model. The region with f 2 lighter than x can n °t be excluded if the LSP is a neutral sparticle 
other than \. One should, though, make sure that the decay of f 2 to the LSP does not destroy 
the predictions of the standard BBN |401 141] . 

We present the restrictions from all the requirements imposed in the M 1( / 2 — tuq plane for 
Aq/Mi/2 = 0, 1, —1, and —2 in Fig. [2j From the relevant data, we observe that the lower 
bound in Eq. (fTU|) is fulfilled for the mass of the CP-odd Higgs boson vtla > 520 GeV and almost 
independently of the other parameters. Note also that, for Aq/Mi/2 = —1 and —2, the bound 
in Eq. (fT2|) is violated for M]/ 2 < 400 GeV and, consequently, does not appear in the relevant 
diagrams. It is obvious that, for all the Aq/ 'M]y 2 's considered in Fig. [21 we are left with no region 
allowed by all the restrictions of Sec. [2j This is due to the fact that the constraint in Eq. (J7J), 
which necessarily holds if x 1S the LSP, is nowhere fulfilled simultaneously with the bound in 
Eq. ([9]). Note, finally, that the upper bound on Mi/ 2 from the lower bound in Eq. (|11|) is also 
nowhere satisfied simultaneously with the bound in Eq. (J9j) for the values taken for Aq/Mi^ in 
Fig. [2j However, this conflict is less serious since, as we already explained in Sec. El the lower 
bound in Eq. (jlip is considered here only as an optional constraint. 

The constraint in Eq. ([7]) is, in general, satisfied in two well-defined distinct regions in the 
diagrams of Fig. [21 which are: 

• The region to the left of the almost vertical part of the line corresponding to the upper 
bound on Mi/ 2 from Eq. ©, where the neutralino annihilation via the s-channel exchange 
of a CP-odd Higgs boson A is by far the dominant (co) annihilation process. However, 
this region is excluded by the constraints in Eqs. ([9]) and (|10p . On the other hand, it is 
well known - see e.g. Refs. [H [16] - that this region is extremely sensitive to variations 
of rrib(Mz)- Indeed, we find that, as rrib(Mz) decreases, the ^4-boson mass rriA increases 
and approaches 2m^. The A-pole neutralino annihilation is then enhanced and f2^/i 2 is 
drastically reduced causing an increase of the upper bound on M.\m. However, even if we 
reduce nib(Mz), we do not find any ^4-pole neutralino annihilation region which is allowed 
by the requirements in Eqs. ([U]) and fjlOj) . 

• The narrow region which lies just above the light gray area with r 2 lighter than the 
neutralino, where bino-stau coannihilations [3l [M] take over leading to a very pronounced 
reduction of Q^h 2 . A large portion of this region survives after the application of all the 
other requirements of Sec.[2]except for that in Eq. Q. To get a better understanding of this 
region, we can replace the parameter mo by the relative mass splitting Af 2 = (nif 2 — m^)/m^ 
between x an d the lightest stau, which controls the strength of bino-stau coannihilations. 
The coannihilation region then approximately corresponds to Af 2 = — 0.25. It is evident 
from Fig. [2] that the slope of the boundary line with Af 2 = increases as Aq/M±/2 moves 
away from zero in both directions. Note that this slope in our model turns out to be larger 
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Figure 2. The restrictions in the M]/ 2 — thq plane for various values of Ao/M±/ 2 indicated in 
the graphs. The conventions adopted are described in Fig. [TJ 



than the one obtained in other versions of the CMSSM - cf. Ref. [3] - with lower values of 
tan/3. As a consequence, small variations of uiq or M^/ 2 lead, in our model, to more drastic 
variations in Af 2 . 

To investigate further whether the incompatibility between the constraints in Eqs. ([7]) and 
([9]) extends to all possible Aq/Mi/ 2 s, we focus on the coannihilation regime and construct the 
regions allowed by all the restrictions of Sec. [2] in the M 1 / 2 — Aq/JV^m plane for Af 2 = 0. We 



depict our results in Fig. The choice Af 2 = ensures the maximal possible reduction of 
VLy>h? due to the x — T2 coannihilation. So, for a given value of Aq/Mi/ 2 , the maximal Mij 2 
or rriy. allowed by Eq. ([7|), which holds under the assumption that the neutralino is the LSP, 
corresponds to Af 2 = 0. We find that, for Aq/M 1 / 2 < 0, processes with f 2 f 2 in the initial 
state and IF ± VF =F , W ± H ZfZ in the final one become more efficient (with a total contribution 
to the effective cross section of about 14 to 22% as Aq/Mi/ 2 decreases from to —2) and so 
coannihilation is strengthened and m^s larger than in the Aq/Mi/ 2 > case are allowed by 
Eq. ©. The overall maximal M 1/2 ~ 1575 GeV or ~ 722 GeV allowed by Eq. is 
encountered at Aq/Mi/ 2 ~ —2 yielding BR(£? S — > fj, + /j,~) = 1.82 x 10 -8 . Comparing the above 
upper bound on M^j 2 with the corresponding one in Fig. 3 of Ref. [22] (represented by a solid and 
two dotted black lines), we observe that here the bound is considerably enhanced in the region 
of low as well as the region of large values of Aq/Mi/ 2 since we do not consider the constraint 
from the lower bound on Sa^ from the r-based calculation. However, it always remains smaller 
than the lower bound on Mi/ 2 derived from Eq. © - note that for 3.7 < Aq/Mi/ 2 < 3.9 the 
overall lower bound on M\j 2 is derived from Eq. ( 1101 ). Indeed, the smallest lower bound on 
M 1/2 = 1306 GeV or 



590 GeV is found at A /M 1/2 ~ 2 yielding Q^h 2 = 0.15. Note that 
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Figure 3. The restrictions in the 
Mi/2 - Aq/Mx/2 plane for A f2 = 
following the conventions of Fig. [TJ 



increasing A^ 2 within the range — 0.25 does not alter the boundaries of the various constraints 
in any essential way, except the solid line which is displaced to the left shrinking, thereby, the 
area allowed by Eq. ([7]) considerably. Needless to say that the more stringent optional upper 
bound on M\/2 from the lower bound in Eq. (jlip is also not compatible with the constraint in 
Eq. ([9]). Consequently, for every Aq, there is no range of parameters simultaneously allowed by 
all the constraints and, therefore, x can De now excluded as a CDM particle in our model. 

The exclusion, in our model, of x as a CDM candidate, resulting from the incompatibility 
between Eqs. ([7]) and (|9|), is further strengthened if one tries to reconcile Eqs. d7J and (fT3|) . 
Indeed, the tension between the neutralino CDM and the new data on is quite generic within 
the CMSSM since the fulfillment of Eq. (|13p requires a very heavy SUSY spectrum, which leads 
to conflict with Eq. ([7]) - cf. Refs. [Ml E7J- It would be interesting to investigate this issue in our 
model, which yields large values of tan f3 and so Eq. (|13p can be possibly satisfied with a lighter 
SUSY spectrum relative to other versions of the CMSSM with lower tan /3's. Our results are 
presented in Fig. HI where we draw rrih (bold lines) versus for Af 2 ~ and Aq/Mi/ 2 = —2, 
— 1, 0, 1, and 3.9. In the table included in Fig. HI we also list the minimal m^-'s, m^| m i n , for 
which the inequality in Eq. © and the lower bound in Eq. (|10p are satisfied for given values of 
A0/M1/2 as well as the corresponding ranges of Q%h 2 as varies from m^lmm to about 1 TeV. 
Since we take Af 2 ~ 0, the derived Q-c/i 2 takes its minimal possible value. In this plot, we also 
depict by a dotted line the value of m^lmm for all possible Aq/Mi^'s and A^ 2 ~ 0. This line 
terminates at Aq/M±/2 — 3.9 since beyond this value the stability of the electroweak vacuum 
fails. The overall minimal ~ 119 GeV is encountered at Aq/Mi/2 — 2 and ~ 590 GeV. It 
is interesting to note that m;, increases with (or Mi/ 2) an d as Aq decreases and, eventually, 
this mass enters inside the gray region in Fig. which is preferred, at la, by the recent LHC 
searches. However, it is obvious from Fig. U] and the values of J7y/t 2 in the table included in this 
figure that Eqs. (J7J) and (|13p are incompatible even in our model. 

Departure from Af 2 = is not expected to alter drastically our predictions as regards the 
value of rrih since m-h depends crucially on M1/2, but only mildly on m,Q. Moreover, one can 
deduce from the slop of the left boundary of the gray regions in Fig. [2] that, increasing Af 2 , 
smaller Mi/ 2 's and, therefore, slightly lighter Higgs masses are permitted by Eq. ©. Note that, 
since x cannot be the LSP, both signs of Af 2 are possible. In particular Af 2 > [Af 2 < 0] 
corresponds to x [^2] being the LOSP. The critical case Af 2 ~ gives the minimal possible relic 
abundance of the LOSP in both cases due to the coannihilation effect and, therefore, the maximal 
possible mass of the LSP if this is produced mainly non-thermally - see Refs. [40 [ 1411 H2] . 

Having in mind mostly the latter possibility, we proceed in the presentation of our predictions 
for the sparticle and the Higgs boson spectrum of our model, which may be observable at the 
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Figure 4. The variation of rrih as a function of for Af 2 ~ and various Aq/Mi^'s indicated 
on the curves and in the table included in this figure. In this table, listed are also the minimal 
m%'s for which the inequality in Eq. ([9]) and the lower bound in Eq. ()10p are satisfied and the 
corresponding ranges of Qj^h? as Tn^ increases from fn^lmin 

to about 1 TeV. The minimal nih 
for all the values of Aq/Mu^ is also depicted by a dotted line. The part of the region of Eq. ([13]) 
preferred by the recent LHC data which lies in the panel is painted light gray. 




LHC. In Table [T] we list the model input and output parameters, the masses in GeV of the 
sparticles - neutralinos x, x%i xl], x% charginos Xn xf? gluinos g, squarks t\, t2, &i, &2 5 &L> ur, 
dz, d,R, and sleptons fi, T2, u r , bl, &r, v e - and the Higgs bosons (h, H, H^, A), and the values 
of the various low energy observables for Af 2 ~ 0, Aq/Mi/ 2 = —3, —2, —1, 0, and 2 and for 
the minimal Mj/2 allowed by Eq. ([9]) in each case. We consider the squarks and sleptons of the 
two first generations as degenerate. From the values of the various observable quantities, we see 
that the bound in Eq. © and the optional lower bound in Eq. (flT]) are violated. So, the lightest 
neutralino cannot be the LSP. It is also very interesting to observe that the predicted values of 
rrih lie close or even inside the range in Eq. ([13]) favored by Ref. [25] - cf. Refs. [54"] IFF]. 

The deviation from YU can be estimated by defining [19] the relative splittings 8h b and 5h T 
at Mqut through the relations: 

_h b -h t 2c _h t -h T 

oh h = — = -— — = -6h T = — (14) 

h t 1 + c ht 

Along the dotted line of Fig. 0] the ranges of the parameters c, Sh T , 5hb, and tan/3 are 

0.148 < c < 0.16, 0.26 < 6h T = -5h b < 0.28, 56.2 < tan/3 < 56.9. 

Let us underline that, although the required deviation from YU is not so small, the restrictions 
from YU are not completely lost since tan (3 remains large - close to 60 - and that the deviation 
from exact YU is generated within well- motivated SUSY GUTs described in Ref. |16j . 

It is worth emphasizing that our results do not invalidate the x candidacy for a CDM 
particle in all versions of the CMSSM with Yukawa quasi-unification. This is because the 
(monoparametric) condition of Eq. ([1]), which we considered here, is only a simplified case of the 
Yukawa quasi-unification conditions shown in Eq. (15) of Ref. [16], which depend on one real and 
two complex parameters. Actually, the investigation of the viability of x as a CDM candidate 
within the CMSSM with Yukawa quasi-unification conditions more complicated than the one 
in Eq. ([1]) derived from the GUT models of Ref. [16] is under consideration. Alternatively, our 
present model may be perfectly consistent with data if we avoid the restriction from Eq. ([7]) by 
assuming that the LSP is the axino |41l |4"2] with mass a little lower than and that the reheat 
temperature is adequately low. 



Table 1. Input and output parameters, masses of the sparticles and Higgs bosons, and values 
of the low energy observables for Af 2 ~ 0, five values of Aq/ 'Mi/ 2 , and the minimal Mi/ 2 . 



Input parameters 



A Q /M 1/2 


-3 


-2 


-1 





2 


c 


0.1589 


0.1592 


0.1585 


0.153 


0.1475 


M 1/2 /GeV 


2355.02 


2019.6 


1744.65 


1530.3 


1317.65 


mo/GeV 


4292.06 


2760.01 


1691.01 


1132.04 


1542.88 


Output parameters 


tan j3 


57.1 


56.9 


56.3 


56.2 


56.2 


100(5 /i t (M G ut) 


27.4 


27.4 


27.3 


26.5 


25.7 


fi/GeV 


3704.64 


2755 


2059 


1588 


1250 


Masses in 


GeV of sparticles and Higgs 


bosons 




X 


1089.2 


926.5 


794.2 


692.5 


595.2 


A2 


2087.2 


1777.7 


1524.2 


1325.4 


1130.1 


X°s 


3688.5 


2747.6 


2057.3 


1588.9 


1254.8 


xl 


3789.7 


2750.3 


2062.7 


1600.5 


1279.9 


xt 


3690.1 


2750.6 


2062.8 


1600.4 


1279.6 


xt 


2087.3 


1777.8 


1524.3 


1325.5 


1130.2 


9 


5190.9 


4454 


5190.9 


3388.6 


2981.2 


ti 


4336.3 


3608.2 


3094.3 


2752.2 


2567.8 


ti 


3593.2 


3084.5 


2709.5 


2449.2 


2303.1 


h 


4514.1 


3653.4 


3097.4 


2747.9 


2575.7 


h 


4310.6 


3561.2 


3005.7 


2644.1 


2508.8 


UL 


6215.9 


4786 


3815.1 


3231.7 


3051.1 


UR 


6047.6 


4624.8 


3661.4 


3090.5 


2941.3 


d L 


6216.3 


4786.5 


3815.8 


3232.5 


3052 


<Ir 


6026.1 


4604.4 


3641.9 


3072.7 


2927.7 


h 


3447.3 


2413.2 


1721.0 


1354.4 


1436.4 


T2 


1089.9 


927.1 


794.4 


692.5 


595.7 


V T 


3443.9 


2407.7 


1712.4 


1343.3 


1430.2 




4582.6 


3085.1 


2070.4 


1544.9 


1789.6 


e-R 


4389.7 


2869.5 


1821.4 


1278.5 


1626.4 




4581.6 


3083.7 


2068.6 


1542.6 


1787.5 


h 


126.17 


124.3 


122.68 


121.15 


119.30 


H 


1463.72 


1334.6 


1181.53 


1012.4 


730.78 




1466.38 


1337.6 


1185.04 


1016.55 


736.91 


A 


1463.99 


1334.9 


1182.00 


1013 


732 


Low energy observables 


10 4 BR(fe -)■ s 7 ) 


3.25 


3.23 


3.22 


3.23 


3.35 


10 8 BR(B S 


1.08 


1.08 


1.08 


1.08 


1.08 




0.929 


0.915 


0.893 


0.856 


0.736 


lO 10 ^ 


0.565 


1.09 


2.04 


3.27 


3.4 


%/i 2 


0.301 


0.219 


0.167 


0.152 


0.151 



4. Conclusions 

We performed a revised scan of the parameter space of the CMSSM with [i > applying a 
suitable Yukawa quasi-unification condition predicted by the SUSY GUT model of Ref. |16j . 
which has been constructed in order to remedy the 6-quark mass problem arising from exact 
YU and universal boundary conditions. We took into account updated constraints from collider 
and cosmological data. These constraints originate from the CDM abundance in the universe, 
B physics (6 — > sj, B s — > and B u — > rz/), 5a ^, and m^. Although the neutralino- 

stau coannihilations drastically reduce the neutralino relic abundance and, thus, enhance the 
upper bound on implied by the assumption that the neutralino is a CDM particle, they do 
not quite succeed to bring it to an acceptable level compatible with the lower bound on my. 
induced by BR (B s — > Therefore, - contrary to our findings in Ref. [22J - x is excluded 

as CDM particle by the combination of the constraints from BR(i? JS — > and CDM. As 

a consequence, the model can become consistent with observations only if the LSP is a SUSY 
particle other than the neutralino. This could be the axino or the gravitino and can account 
for the present CDM abundance in the universe. It is interesting to note that, in this case, the 
lowest predicted is enhanced and gets closer to the range favored by the recent preliminarily 
results announced by LHC. 
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